o . 
o ■ 
o ■ 
<N : 

^3 : 

'. Low-dimensional BEC 

O ■ 

m ; 

F.J. Sevillai, M. Grether^, M. Fortes^ M. de Llano^ O. Rojo'^, 

' pH'l M.A. Soli's^ and A. A. Valladares'^ 

o ■ 

y [ ^Instituto de Fisica, UN AM, Apdo. Postal 20-364, 01000 Mexico, DF, Mexico 

■ '^Facultad de Ciencias, UNAM, 04510 Mexico, DF, Mexico 
'^Instituto de Investigaciones en Materiales, UNAM, 04510 Mexico, DF, Mexico 

^ ■ ^PESTIC, Secretana Academica & CINVESTAV, IPN, 04430 Mexico, DF, Mexico 
^ . 

a: 

c ■ 

Q ' The Bose-Einstein condensation (BEC) temperature Tc of Cooper pairs 

, (CPs) created from a general interfermion interaction is determined for a 

\ linear, as well as the usually assumed quadratic, energy vs center- of -mass 

^ ■ momentum dispersion relation. This explicit Tc is then compared with a 

. widely applied implicit one of Wen & Kan (1988) in d = 2 + e dimensions, 

lO \ for small e, for a geometry of an infinite stack of parallel (e.g., copper- 

' oxygen) planes as in, say, a cuprate superconductor, and with a new result 

■ for linear- dispersion CPs. The implicit formula gives Tc values only slightly 
O : lo.er tkan tkose of ike e.pi.c, formula for iy„cal cu.raie ^r.ra.U„. 

■ PACS numbers: 74.20.Fg; 64.90+b; 05.30.Fk; 05.30.Jp 

^ ^ \ Bose-Einstein condensation (BEC) of Cooper pairs (CPs) can lead to 

j-H ■ a phase transition (even in 2D) in any many-fermion system dynamically 

O ■ capable of forming CPs. This transition could be the origin of "exotic" 

?^ \ superconductivity in the quasi-2D cuprates and in the quasi- ID organo- 

' metallic (Bechgaard) salts, as well as of the superfluidity in liquid ^He or in 

^ . trapped Fermi gases in 3D. 

^ \ The familiar BEC formula for the transition temperature is 

Tc ~ -i.Zlh^nl^ ImBkB, (1) 

with ub the number density of bosons of mass and the Boltzmann 
constant. This is a special case of the more general expressionS valid for any 
space dimensionality c? > and any boson dispersion relation ek = Cg 
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with s > and Cs constant, given by the explicit Tc-formula 

s/d 



sT{d/2) {2TTynB 
27r'^/2r(d/s)5d/.(l) 



(2) 



If fi{T) is the boson chemical potential and eA»(^)/'=sT ^ ^ ^j^g fugacity, 
ga{z) = Y^'iZi z'' / 1"^ are the Bose integrals. For z = 1 and a > 1 the ga{^) 
is just C('^)) the Riemann Zeta-function of order a which is finite for a > 1 
and infinite for cr = 1, while the series ga{^) diverges for all o" < 1. For 
s = 2, C2 = h^/2mB, and since ({3/2) ~ 2.612, this leads to the usual BEC 
Tc-formula (||). Since 5^/2(1) diverges for all (i/2 < 1, Tc = for all d < 2. 
This follows from the boson number equation 

AT = No{T) + J2 [e^^^-'^(^)>/'=^^ - 1] (3) 



c 



where Nq{T) is the number of bosons in the K = state. At T = 
both Nq(Tc) and the boson chemical potential fJ-{Tc) virtually vanish so that 
replacing 

^ ^ (L/2^)'^/ d'^fc = (L/2vr)'^— — / dkk'^^ (4) 
^ Jo+ r((i/2) Jo+ 

in (^) eventually yields (|2|) where ub = N/L'^. 

The fact that a CP can have a linear (s = 1), as opposed to the usual 
quadratic fs = 2), dispersion relation was mentioned as far back as 1964 by 
Schrieffer.i p. 33, for the BCS model interaction in 3D. This was recently 
confirmedB to be the case for both 2D and 3D under a very general inter- 
fermion interaction for any coupling provided the fermion number-density 
is nonzero, i.e., in the presence of a Fermi sea. The CP dispersion relation 
becomes quadratic only in the extremely dilute (or vacuum) limit where the 
CPs are just the so-called "local pairs." For any sizeable fermion density the 
nonnegative CP excitation energy ek = Aq — Ax behaves like ~ a{d)hvFK , 
where Ax (not to be confused with the BCS gap A) is the (positive) binding 
energy of a CP of center-of-mass momentum (CMM) hK, vp = fikp/m and 
kp the Fermi velocity and wavenumber, respectively, m the fermion effective 
mass, while a{d) = 2/7r and 1/2 in 2D and 3D, respectively, precisely as 
established^ previously for the BCS model interaction. For linear dispersion 
s = 1, Ci = a{d)hvF, Tc = from (§) for all d < 1 only — and Tc > for all 
d > 1, which is precisely the range of dimensionalities for all known super- 
conductors if one includes the quasi-lD organo-metallic Bechgaard salts.! 
Using the interpolation a{d) = (7/2 - 6/tt) + (S/vr - 13/4)(i+ (3/4 - 2/7r)d2, 
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which correctly reduces to 1, 2/7r and 1/2 in ID, 2D and 3D, respectively. 
Fig. 1 graphs (|2|) for s = 1 and 2 (in units of the Fermi temperature Tp) 
vs d — if one imagines all the fermions in the initial, interactionless many- 
fermion system paired into CPs of mass ms = 2m. The particle number 
density of the original fermions is n = kj,/2'^-^TT'^/^d T{d/2) and equals 2nB, 
and we have used tv^kp/lm = ^mvp = Ep = knTp. These curves are upper 
hounds to the Tc from a more realistic modelQ where chemical equilibrium 
allows only a fraction of all fermions to be actually bound into pairs. 
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Fig. 1. Left: dimensionality, d, dependence of the critical BEC transition 
temperature Tc according to explicit (^) (in units of the Fermi temperature 
Tp) as explained in text. Lower (full) curve is for s = 2; upper (dashed) 
curve for s = 1. Shaded areas refer to empirical data from Ref. [7]. Right: 
dots refer to results from implicit (^) below, as explained just after (^). 

A rather general interfermion interaction is the S'-wave attractive sepa- 
rable potential whose double Fourier transform is 

Vpq = -{vo/L^)9p9q- (5) 

Here L is the size of the "box" confining the many-fermion system, vq > 
is the interaction strength and given, e.g.,i by (1 +p^/Po)~^^^ where po 
is the inverse range of the potential. Hence, po ^ oo implies = 1 and 
corresponds to the attractive contact (or delta) potential V{r) = —vo5{r), 
while Po = kp implies a range of order of the average interfermion spacing, 
etc. If Qp = 9{h.u)D + f^F — p^ /2m), with 6{x) the unit step function, ^ 
becomes the BCS model interaction where a;^) is the Debye frequency and 
Hp the fermionic chemical potential that becomes Ep for T = = uq. 

Using a renormalized CP equatior^ whose coupling depends only on the 
two-body binding energy B2P the CP excitation energy ek = ^0 — '^K for 
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zero range was obtained numericallyil as an exact curve that for very small 
B^jEp is virtually linear, i.e., ek IfivpKjii . It is only in the dilute 
limit {vp or Ep 0) that £k tends asymptotically to the exact quadratic 
1)?K'^ /2{2m) for any coupling. Assuming ns = n/2 and tjib = 2m to 
introduce the temperature scale Tp as before. Fig. 2 shows the BEC Tc's of 
a pure gas of unbreakable CPs. Significantly, Tc is no longer zero in 2D — as 
would be predicted in a BEC picture by a quadratic relation appropriate for 
"local-pair" CPs in vacuum, a result that wrongly suggests that BEC cannot 
apply for quasi-2D cuprate superconductors. 

More accurate BEC Tc's should include refinements such as non-5'-wave 
interactions, allowing for unpaired fermions in a more realistic binary boson- 
fermion mixture modelQ — and most importantly, CP-fermion interactions 
that link0'0 the BEC condensate fraction temperature-dependence with 
that of the BCS fermionic energy gap, among other corrections. 

The linear dispersion relation of a CP should not be confused with the 
linear dispersion of Anderson-Bogoliubov-Higgs (ABH) many-body excita- 
tion phonon-like modes. Collective modes in a superconductor were studied 
since the late 1950's bvseveral workers. A more recent treatment for ID, 
2D and 3D is availablell3 which confirms the linear ABH form hvpKj^fd for 
d =1, 2 or 3 in the zero-coupling limit. Our CPs are taken as "bosonic" 
even though they do not obey (Ref. [2] p. 38) Bose commutation relations. 
This is because for a given K they have indefinite occupation number since 
for fixed K there are (after the thermodynamic limit) an indefinitely large 
number of allowed (relative wavenumber) k values, so that — for any cou- 
pling and thus any degree of overlap between them — CPs do in fact obey 
the Bose-Einstein distribution from which BEC is determined. By contrast, 
ABH phonons (like photons or plasmons, etc.) cannot suffer a BEC as their 
number is always indefinite. The number of CPs, on the other hand, is fixed 
at half the number of (pairable) fermions if all of these are imagined paired 
at a given temperature and coupling. 

To model cuprate superconductors consider the bosons confined to an 
infinite set of planes stacked along the z-direction, parallel to each other with 
equal spacing c between adjacent planes. The BEC transition temperature 
formula, for s = 2 bosons in each plane, is the implicit Tc-formulaE^ 

k T = ^imB-zph^c 

" mnHMBC^kBTMQ/h^] ' ^ ^ 

and is valid in 2 + e dimensions where e is small and given by 

e^2[ln{nB-3DMBC^/mB)]-\ (7) 

This expressly vanishes as Mbc^ —>■ oo, as it should. Here ub-zd = N/L^ 
while u{t) = 1 - t + 0{t^) ~ 1 if t << 1, where t = H^/MBC^kBT = 
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Fig. 2. 2D critical BEC temperatures Tc (in units of Tp) for five coupling 
values B2/EP for a pure boson gas of unbreakable CPs, determined by the 
Nq{Tc) = = /u(rc) solution of (^) inserting exact numerical CP dispersion 
curves, using (^). The value 0.702 marked with a square corresponds to the 
Tc/Tp value at zero coupling (see also Ref. [1]) and contrasts with the well- 
known result Tc/Tp = in 2D for infinite coupling where £k = h'^K'^ /2{2m) 
exactly. Shaded areas as in Fig. 1. 

n^/2mc^{MB/mB)kBT. Using h^/mkB = 88,419 K with m the 
electron mass, and c = 12 A, this inequality is well satisfied for the 
higher cuprate transition temperatures, today ranging up to 164 K, since 
Tc >> 307K/(iW^/m^) as typically MB/rnB can range from 10^ to 10^. 
Clearly, for Mbc — > (infinitely separated planes and/or perfect confine- 
ment to the z-direction in each plane) Tc vanishes as it should in 2D. As 
state, this Tc-equation is implicit or transcendental, unlike the simpler ex- 
plicit Tc equations (J) and and has been used for varied purposes by 
numerous authors0B0 — though only for quadratic dispersion bosons. 
We have generalized @ for any s > and found 



kBTc = Cs [27rnB-3D5c/r(2/s)<72/.(e-^'/^^^^''=^^= 



s/2 



(8) 



Thus, an exact (again, implicit) equation for Tc is obtained for any s > 0. 
For s = 2 and C2 = 'h? /2mB we recover (^. In Fig. 1 (right) we plot results 
for both values of s as points, for Mb/ttib = 10^ and 10^. They can be seen 
to differ very slightly from the results for s = 2 or s = 1 bosons in d = 2 -|- e 
dimensions, with e small, that came directly from (^) which, moreover, is 
valid for all d > 0. 

To compare our results, consider the Berezinskii-Kosterlitz-ThoulesstS 
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transition temperature formula 

fcBTf = 5^ (9) 
2 niB 

valid in 2D, and assume as before that hb = n/2 = kp/An and niB = 2m. 
This gives T^^^/Tp = 1/8 = 0.125 and is displayed as a square in Fig. 1. 

In conclusion, BEC T^s related to a pure gas of unbreakable composite 
linear-dispersion bosons were calculated in d = 2 for Cooper pairs formed 
via a general separable potential and whose coupling for any CMM is char- 
acterized solely by its two-body binding energy in vacuum. A Tc-formula 
valid for any dispersion relation of the form ex = Cg with s > and Cs 
constant is deduced for the BEC of a pure gas of unbreakable bosons con- 
fined to move in infinitely-many identical planes parallel to each other. It is 
a peculiar implicit Tc-equation valid in 2 + e where e is small and vanishes, 
as it must, when the inter-plane spacing or the boson mass in the direction 
perpendicular to the planes diverges. However, for either s = 2 or 1 we 
find results that are just slightly different than those of the explicit BEC 
Tc-formula which is valid more generally for any d > 0. 
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CONACyT (Mexico) # 27828 E, DGES (Spain) # PB95-0492 and FAPESP 
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The Bose-Einstein condensation (BEC) temperature Tc of Cooper pairs 
(CPs) created from a general interfermion interaction is determined for a 
linear, as well as the usually assumed quadratic, energy vs center-of-mass 
momentum dispersion relation. This explicit Tc is then compared with a 
widely applied implicit one of Wen & Kan (1988) in d = 2 -\- e dimensions, 
for small e, for a geometry of an infinite stack of parallel (e.g., copper- 
oxygen) planes as in, say, a cuprate superconductor, and with a new result 
for linear-dispersion CPs. The implicit formula gives Tc values only slightly 
lower than those of the explicit formula for typical cuprate parameters. 

PACS numbers: 74.20.Fg; 64.90+b; 05.30.Fk; 05.30.Jp 

Bose-Einstein condensation (BEC) of Cooper pairs (CPs) can lead to 
a phase transition (even in 2D) in any many-fermion system dynamically 
capable of forming CPs. This transition could be the origin of "exotic" 
superconductivity in the quasi-2D cuprates and in the quasi- ID organo- 
metallic (Bechgaard) salts, as well as of the superfluidity in liquid ^He or in 
trapped Fermi gases in 3D. 

The familiar BEC formula for the transition temperature is 

Tc - ^.2.in'^njl^ /ruBkB, (1) 

with ub the number density of bosons of mass rriB and kB the Boltzmann 
constant. This is a special case of the more general expression-^ valid for any 
space dimensionality d > and any boson dispersion relation ex = Cg K" 
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with s > and Cg constant, given by the explicit Tc-formula 

1 



sT{d/2) {2TTYnB 
27rrf/2r(d/.)ff,/,(l) 



(2) 



If /i(T) is the boson chemical potential and eM'^^^^^'^ = z the fugacity, 
ga{z) = Y^fZi ^ jV' are the Bose integrals. For z = \ and c > 1 the ga^S) 
is just C(c)j the Riemann Zeta-function of order u which is finite for cr > 1 
and infinite for cr = 1, while the series ^^(l) diverges for all cr < 1. For 
s = 2, C2 = n'^/2mB, and since C(3/2) ~ 2.612, this leads to the usual BEG 
Tc-formula (1). Since 5^/2(1) diverges for all d/2 < 1, = for all d<2. 
This follows from the boson number equation 



c 



where Nq{T) is the number of bosons in the K = Q state. At T = 
both No{Tc) and the boson chemical potential fJ,{Tc) virtually vanish so that 
replacing 



y — > {L/2Trf d'^k= (L/2Tr)'^—— dkk'^-^ (4) 



in (3) eventually yields (2) where ub = N/L'^. 

The fact that a CP can have a linear {s = 1), as opposed to the usual 
quadratic (s = 2), dispersion relation was mentioned as far back as 1964 by 
Schrieffer,^ p. 33, for the BCS model interaction in 3D. This was recently 
confirmed^ to be the case for both 2D and 3D under a very general inter- 
fermion interaction for any coupling provided the fermion number-density 
is nonzero, i.e., in the presence of a Fermi sea. The CP dispersion relation 
becomes quadratic only in the extremely dilute (or vacuum) limit where the 
CPs are just the so-called "local pairs." For any sizeable fermion density the 
nonnegative CP excitation energy ck = Aq — Ak behaves like a{d)hvpK , 
where Ak (not to be confused with the BCS gap A) is the (positive) binding 
energy of a CP of center-of-mass momentum (CMM) TiK, vp = Tikp/m and 
kp the Fermi velocity and wavenumber, respectively, m the fermion effective 
mass, while a{d) = 2/tt and 1/2 in 2D and 3D, respectively, precisely as 
established^ previously for the BCS model interaction. For linear dispersion 
s = 1, Ci = a{d)hvp, Tc = from (2) for all d < 1 only — and Tc > for all 
c? > 1, which is precisely the range of dimensionalities for all known super- 
conductors if one includes the quasi-lD organo-metallic Bechgaard salts. ^ 
Using the interpolation a{d) = (7/2 - 6/7r) + (S/tt - 13/4)d+ (3/4 - 2/7r)d^ 
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which correctly reduces to 1, 2/7r and 1/2 in ID, 2D and 3D, respectively. 
Fig. 1 graphs (2) for s = 1 and 2 (in units of the Fermi temperature Tp) 
vs d — if one imagines all the fermions in the initial, interactionless many- 
fermion system paired into CPs of mass tub = 2m. The particle number 
density of the original fermions is n = kp/2'^^'^'7r'^/'^d T{d/2) and equals 2n5, 
and we have used h^k^/2m = ^mvp = Ep = UbTp. These curves are upper 
hounds to the Tc from a more realistic model^ where chemical equilibrium 
allows only a fraction of all fermions to be actually bound into pairs. 




Fig. 1. Left: dimensionality, d, dependence of the critical BEC transition 
temperature Tc according to explicit (2) (in units of the Fermi temperature 
Tp) as explained in text. Lower (full) curve is for s = 2; upper (dashed) 
curve for s = 1. Shaded areas refer to empirical data from Ref. [7]. Right: 
dots refer to results from implicit (6) below, as explained just after (7). 

A rather general interfermion interaction is the 5- wave attractive sepa- 
rable potential whose double Fourier transform is 

Vpq = -{vo/L^)gpgq- (5) 

Here L is the size of the "box" confining the many-fermion system, > 
is the interaction strength and given, e.g.,^ by (1 + p^ / Pq)^^ ^"^ where pQ 
is the inverse range of the potential. Hence, pQ ^ oo implies = 1 and 
corresponds to the attractive contact (or delta) potential V{r) = —VQ6{r), 
while po = kp implies a range of order of the average interfermion spacing, 
etc. If Qp = 6{Tiujd + lip — p^ /2m), with 6{;x) the unit step function, (5) 
becomes the BCS model interaction where cjj) is the Debye frequency and 
^ip the fermionic chemical potential that becomes Ep for T = = vq. 

Using a renormalized CP equation^ whose coupling depends only on the 
two-body binding energy B2,^ the CP excitation energy £k = Aq — Ak for 
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zero range was obtained numerically^ as an exact curve that for very small 
B2/EP is virtually linear, i.e., ek 2TivpK /tt. It is only in the dilute 
limit {vp or Ep 0) that ex tends asymptotically to the exact quadratic 
h?K^ /2{2m) for any coupling. Assuming hb = n/2 and niB = 2m to 
introduce the temperature scale Tp as before, Fig. 2 shows the BEC Tc's of 
a pure gas of unbreakable CPs. Significantly, Tc is no longer zero in 2D — as 
would be predicted in a BEC picture by a quadratic relation appropriate for 
"local-pair" CPs in vacuum, a result that wrongly suggests that BEC cannot 
apply for quasi-2D cuprate superconductors. 

More accurate BEC Tc's should include refinements such as non-iS-wave 
interactions, allowing for unpaired fermions in a more realistic binary boson- 
fermion mixture model^ — and most importantly, CP-fermion interactions 
that link^'^'^^ the BEC condensate fraction temperature-dependence with 
that of the BCS fermionic energy gap, among other corrections. 

The linear dispersion relation of a CP should not be confused with the 
linear dispersion of Anderson-Bogoliubov-Higgs (ABH) many-body excita- 
tion phonon-like modes. Collective modes in a superconductor were studied 
since the late 1950's by several workers. A more recent treatment for ID, 
2D and 3D is available^^ which confirms the linear ABH form TivpK /^/d for 
d =1, 2 or 3 in the zero-coupling limit. Our CPs are taken as "bosonic" 
even though they do not obey (Ref. [2] p. 38) Bose commutation relations. 
This is because for a given K they have indefinite occupation number since 
for fixed K there are (after the thermodynamic limit) an indefinitely large 
number of allowed (relative wavenumber) k values, so that — for any cou- 
pling and thus any degree of overlap between them — CPs do in fact obey 
the Bose-Einstein distribution from which BEC is determined. By contrast, 
ABH phonons (like photons or plasmons, etc.) cannot suffer a BEC as their 
number is always indefinite. The number of CPs, on the other hand, is fixed 
at half the number of (pairable) fermions if all of these are imagined paired 
at a given temperature and coupling. 

To model cuprate superconductors consider the bosons confined to an 
infinite set of planes stacked along the 2;-direction, parallel to each other with 
equal spacing c between adjacent planes. The BEC transition temperature 
formula, for s = 2 bosons in each plane, is the implicit Tc-formula^^ 

kBT. = — ^""r^^^''^ — ^, (6) 

mBHMBc''kBT,iy{t,)/hy 
and is valid in 2 + e dimensions where e is small and given by 

e - 2[\ii{nB-^DMBC^/mB)]-^. (7) 

This expressly vanishes as Mbc^ 00, as it should. Here ub-sd = N/L^ 
while u{t) = I - t + 0{t'^) ~ 1 if t « 1, where t = h^/MBC^kBT = 
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Fig. 2. 2D critical BEC temperatures Tc (in units of Tp) for five coupling 
values BijEp for a pure boson gas of unbreakable CPs, determined by the 
Nq{Tc) = = l-i{T(^ solution of (3) inserting exact numerical CP dispersion 
curves, using (4). The value 0.702 marked with a square corresponds to the 
Tc/Tp value at zero coupling (see also Ref. [1]) and contrasts with the well- 
known result Tc/Tp = in 2D for infinite coupling where ex = Tl^K"^ /2{2m) 
exactly. Shaded areas as in Fig. 1. 

/2mc^{MB/mB)kBT. Using h^/mkB = 88,419 K with m the 
electron mass, and c = 12 A, this inequality is well satisfied for the 
higher cuprate transition temperatures, today ranging up to 164 K, since 
Tc >> 307K. / [M B / m b) as typically Mb/itib can range from 10^ to 10^ 
Clearly, for Mbc (infinitely separated planes and/or perfect confine- 
ment to the z-direction in each plane) Tc vanishes as it should in 2D. As 
state, this Tc-equation is implicit or transcendental, unlike the simpler ex- 
plicit Tc equations (1) and (2), and has been used for varied purposes by 
numerous authors"^'^' — though only for quadratic dispersion bosons. 
We have generalized (6) for any s > and found 



kBTc = Cs \27rnB-3DSc/r{2/s)gy,{e-'''/^^'^'^^^'^ 



s/2 



(8) 



Thus, an exact (again, implicit) equation for Tc is obtained for any s > 0. 
For s = 2 and C2 = h'^/2mB we recover (6). In Fig. 1 (right) we plot results 
for both values of s as points, for Mb /tub = 10^ and 10^. They can be seen 
to differ very slightly from the results for s = 2 or s = 1 bosons in d = 2 + e 
dimensions, with e small, that came directly from (2) which, moreover, is 
valid for all c? > 0. 

To compare our results, consider the Berezinskii-Kosterlitz-Thouless^^ 
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transition temperature formula 

^bT^ = ^— — (9) 

valid in 2D, and assume as before that ub = n/2 = kp/Air and niB = 2m. 
This gives T^^^ /Tp = 1/8 = 0.125 and is displayed as a square in Fig. 1. 



In conclusion, BEC Tc's related to a pure gas of unbreakable composite 
linear- dispersion bosons were calculated in c? = 2 for Cooper pairs formed 
via a general separable potential and whose coupling for any CMM is char- 
acterized solely by its two-body binding energy in vacuum. A Tc-formula 
valid for any dispersion relation of the form ex = Cg with s > and Cg 
constant is deduced for the BEC of a pure gas of unbreakable bosons con- 
fined to move in infinitely-many identical planes parallel to each other. It is 
a peculiar implicit Tc-equation valid in 2 + e where e is small and vanishes, 
as it must, when the inter-plane spacing or the boson mass in the direction 
perpendicular to the planes diverges. However, for either s = 2 or 1 we 
find results that are just slightly different than those of the explicit BEC 
Tc-formula which is valid more generally for any d > 0. 
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